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Abstract 

We propose an ansatz which solves the Dyson-Schwinger equation for the real scalar fields in 
Poincare patch of de Sitter space in the IR limit. The Dyson-Schwinger equation for this ansatz 
reduces to the kinetic equation, if one considers scalar fields from the principle series. Solving 
the latter equation we show that under the adiabatic switching on and then off the coupling 
constant the Bunch-Davies vacuum relaxes to the state with the flat Gibbons-Hawking density 
of out-Jost harmonics on top of the corresponding de Sitter invariant out-vacuum. 



I. INTRODUCTION 

The goal of the present paper is to understand the impact of large IR loop corrections on the 
vacuum states in fleld theory on Poincare patch (PP) of de Sitter (dS) space. In [l| one-loop 
correction to the scalar fleld Whiteman function was calculated in PP over the Bunch-Davies 

n 

(BD) state [2]. There are large IR loop contributions even for the very massive flelds. 

They reveal themselves through the particle creation — via the vacuum averages (a^ a) 
and {a a), where a and are annihilation and creation operators. E.g. for the real massive 
scalar fleld theory with the self-interaction one obtains that {a^ ap) oc \og{pri) and 
{flp a-p) oc \og{pri) as the conformal time approaches the future inflnity, — )■ 0. Here p is 
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the modulus of the spatial momentum. Note that in the expanding universe the future infinity, 
7/ — )■ 0, corresponds to the IR limit of the physical momentum, pr] ^ 0. 

Similar IR contributions do appear in other field theories in PP independently of the spin of 
ihe fields and self-interaction potentials, as long as they do not respect conformal invariance 

In |8| t he observations of [ll] were generalized to the other dS invariant states (so called 
a- vacua j9|,[l0|) and to the states containing finite densities of particles. Furthermore, in js] 
kinetic equation was derived. Its solution sums up the leading IR contributions in all loops. 
One of the goals of the present paper is to show explicitly the latter statement, i.e. to derive 
that kinetic equation directly from the Dyson-Schwinger (DS) equation. 

Tree-level Whiteman function for any a-vacuum respects whole dS isometry group 9(], 10| 
even if one restricts field theory to the PP, which covers only half of dS. There are generators 
of dS isometry group which deform PP. As the result, this symmetry is naively broken in the 
vertexes of the loop integrals to a subgroup, respecting only PP. 

However, one can prove* that for the BD state the variation of the loop contributions, under 
those isometry transformations which deform PP, does vanish. Hence, exact Whiteman function 
over the BD state depends only on the dS invariant distance between its two arguments. But 
for the other a-vacua the isometry is broken in loop integrals down to the subgroup in question. 

In curved space-times (or in flat space curvilinear coordinates) various coordinate systems 
frequently cover only their parts. Hence, to do the calculations in such coordinates one has 
to specify suitable conditions at the boundaries of the corresponding patches. Obviously large 
IR effects are sensitive to the boundary conditions. Hence, if one does not perform a careful 
study of the matching between the boundary conditions, he obtains different physical results 
by doing calculations in different coordinate systems. 

In particular it happens that loop contributions in the global dS space are not just large, but 
they are explicitly IR divergent even for the massive fields^ isl, l|. In this respect dS space is 
similar to the QED in strong background electric fields [14]. The presence of such divergences 
shows that the moment when the interactions or background field are switched on can not be 

n 

taken to the past infinity [1]. This puts an abstraction for the dS isometry invariance of the 
correlation functions in global dS and favors the conclusion that cosmological constant should 



* We would like to thank A.Polyakov for telling us the idea of this proof. Some elements of the proof can be 
found in [lll,[12|. 
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be secularly screened by large IR effects. At least with the appropriate choice of the boundary 
conditions, i.e. with those boundary conditions which do not put dS space on the "life support" 

fl. 

The fact that dS isometry is respected in the loops for BD state is a good sign that cosmo- 
logical constant can not be secularly screened in PP. At least if the initial conditions are just 
mild excitations over the a-vacua. But the presence of the large IR effects means that the BD 
state itself gets modified. 

Indeed, the one-loop correction to the BD Whiteman propagator is G^{Z) oc 

log(Z) G^{Z), when the hyperbolic distance is taken to infinity, Z — )■ oo [Ij. This is just the 
Fourier transform of log{pri) corrections. Thus, the factor A^ log{Z) can be big and the loop 
corrections are not suppressed even if A^ is small. 

Thus, the question is what is the dressed state? We show below that this question is related 
to the following one: What is the fate of small density perturbations over a-vacua in the future 
infinity? To address these questions we derive and solve the kinetic equation which describes 
the dynamics of such density perturbations and, as we have mentioned, by product sums the 
leading IR contributions. 

From the solution we see that if one sets BD state as the initial one at past infinity, where 
it is the ground state of the time-dependent free Hamiltonian, this state gets modified even if 
one switches off the coupling constant at future infinity. It will appear that the result of the 
summation of all loops will contain modifications of the BD propagator, which do not vanish 
as A — )■ in the future infinity, but which can not be seen in the free, A = 0, theory. This 
makes dS space quite different from Minkowski or Anti-de-Sitter spaces 151], where adiabatic 
variations of the self-interactions do not change the true vacuum state. 

The phenomenon we observe is a more complicated variant of the following one. Consider 
simple linear oscillator. In the perfectly linear case the oscillator will remain in an excited state 
forever if it was originally in such a state. However, if one will switch on an interaction of the 
oscillator to an external field and then switch it off, the oscillator will relax to the ground state. 
That will happen independently of the type of the interaction or on the type of the external 
field. At the same time of cause the exact state of the system in question, oscillator plus the 
external field, does depend on the type of the interaction. 

In the second section we propose the dS invariant Kadanoff-Baym equation which may be 
suitable to sum the dS invariant IR corrections exactly over the BD state. This section just 
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gives the idea what kind of problem has to be solved. In the third section we derive the kinetic 
equation which does not respect dS isometry, but, unlike full DS equation, is suitable for the 
separation of the IR renormalization form the UV one. The same equation was derived in jsj. 
It was shown there that its collision integral is annihilated by the Gibbons-Hawking density of 
out-Jost states on top of the out-vacuum. The same state annihilates the collision integral of 
the Kadanoff-Baym equation of the second section up to sub leading terms in IR limit. 

To make the paper self-contained we present the general discussion of the scalar fields in 
PP in the Appendix. All the notations, which are not defined in the main text, can be found 
in the Appendix. 



II. TOWARDS INVARIANT KADANOFF BAYM EQUATION FOR BD STATE 

In this paper we are going to study the following field theory: 



Dots here stand for the higher self-interaction terms, which make the theory stable. The reason 
why we are going to consider below formulas only due to the unstable cubic part of the potential 
is just to simplify them. This instability does not affect our conclusions (sl. 

For the BD state the one loop correction to the Whiteman function G |- was calculated in 

|l| (see as well [lol ) . The result for the sum of the tree-level and one-loop contributions in the 
IR limit, Z —7- oo, is as follows: 



^—d^(t)d,(t)+—(t)^ + -^ 



+ ... 



1 



A2 (1 



4 11 



dxx ^ 



-in T 2 



h\x) 



\og{Z) 



G\iZ). 



(2) 



All notations in this formula and in the formulas that follow are given in the Appendix. 

For large enough D the theory in question becomes non-renormalizable. But in the IR 
limit we do not care about UV divergences and renormalizability of the theory in question. 
We assume that all couplings in all equations below take their physical values, i.e. all UV 
divergent (~ log A) or finite (A^) contributions are absorbed into their renormalization. For 
the propagators which have proper Hadamar behavior the UV divergences in dS space are the 
same as in flat one. Below we are keeping track only of the leading large IR contributions. 
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As seen from ([2]), loops are not suppressed in comparison with the tree-level contribution 
for large enough Z . One has to understand what is the result of the summation of the leading 
IR contributions at all loops. The answer on this question can be obtained from the solution 
of the Dyson-Schwinger (DS) equation: 



G{Zxy) = G°(^xy) + X' j [dW] j [dU]G^{Zxw) HZwu) G{Zuy), (3) 

where G{Z) is the matrix of the exact propagators, while is the matrix of the tree- 

level ones. All propagators in ([3]) are the functions of the invariant distance, because we are 
quantizing over the BD state. 

Having in mind the physical and mathematical origin of the large IR effects 8| we have 
simplified the complete system of DS equations in ([3]). We have assumed that the vertex A does 
not receive any new large IR contributions on top of those which are caused by the contributions 
contained in the two-point functions. 

Eq. (I3l) is not suitable for the summation of only large IR contributions A^log(Z), because 
it does not separate the UV from the IR renormalization. One needs an equation which sums 
up only the leading IR contributions and does not even see the contributions which are either 
suppressed by the higher powers of A or even UV divergent (~ A^logA). The proper equation 
is the kinetic one of the next section. However, it does not respect dS isometry, while one would 
like to sum the dS invariant contributions for the BD state. 

One possible variant is as follows. We apply the Klein-Gordon operator to both sides of (j3]) 
to get rid of its dependence on the initial value of the propagator. This operator, when acting 
on the function of Z, is equivalent to + = (Z^ — l)d^ + D Z dz + Tn?- That converts the 
DS equation into an integrodifferential equation of the Kadanoff-Baym form. Recalling that 
the tree-level Whiteman functions, and G^^., solve the homogeneous equation, while the 
Feynman propagators, G\^ and G^__^ solve the inhomogeneous one, we obtain the following 
equation for the Whiteman function G-j^{Zxy)'- 



[^iyC + D ZxY dz^y + m^] G[Zxy] = 
A^ / [dW] G'^[Zxw + «e] G[Zwy + ie sgn{ri^ - rjy)] + 



+A2 j [dW] G\Zxw + sgn{i]^ - G[Zwy 



it\ 



(4) 
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in the limit Zxy oo. However, we do not see that this equation sums up only the leading IR 
terms and nothing else. Possible way to move further is to apply the ansatz G{Z) = f{Z) G^{Z) 
for Z — )■ oo, where f{Z) is slow in comparison with G^{Z). But instead we are going to find 
the stationary IR solution of this equation by taking a different root. 

As the side remark let us mention that it was argued in 1^ that the result of the summation 
of the IR contributions should be the propagator build with the use of the exact Hartle-Hawking 
state. The one which is obtained via analytical continuation from the sphere and constructed 
with the use of the exact Hamiltonian. Obviously such a state depends on the coupling constant 
A. 

The exact state in dS should as well depend on the A, but besides that we encounter a new 
phenomenon, which can not be grasped through the analytical continuation from the sphere. 
We are going to show that the IR stationary solution of (jl]), the one which annihilates its RHS 
up to subleading terms, does not depend on the coupling constant. I.e. even if we start from the 
BD state (ground state of the free Hamiltonian on the sphere) and then adiabatically switch on 
interactions and eventually switch them off the theory relaxes to another state independently 
from the selfinteractions. 

In the next section we will propose the result of the IR dressing of the BD propagator, 
which, however, will not allow us to fix the function f{Z). Because we will be able to find the 
propagator at the stationary state, which is reached as Z — )■ oo, but not how it approaches the 
stationarity in a dS invariant way. We will find the form of its approach to the stationarity 
only in the circumstances when the dS isometry is broken. 

The hint for the stationary propagator comes from the following observations. First, the 
dressed propagator should respect dS isometry. Second, it should annihilate the RHS of (jl]) up 
to the suppressed terms. 



III. SOLUTION OF THE DYSON SCHWINGER EQUATION IN IR LIMIT 

Let us consider small density perturbation over any a-vacuum. Then the dS invariance of 
the propagators is broken even at tree-level, but we still have large IR contributions. To sum 
them up one as well has to solve the DS equation, but this time it does not respect dS isometry. 



Such a DS equation was solved in IR limit in 
a different perspective. 



. In this note we would like to address it from 
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Due to the relation G'^„ + G^, = G"| , + it is convenient to perform the Keldysh 



-+ 

rotation ^ ' ' - -- 



ITj, [18| to the new basis: Dl^{X,Y) = -\ {G\_{X,Y) ^ G\{XX% D^{X,Y) = 
e{r]y - T],) [GO Y) - Gl_iX, Y)] and D^iX, Y) = e{r^, - r]y) [Gl_iX, Y) - +(X, Y)] . 
Here D^'^ are retarded and advanced Green functions. They carry information about the 
quasi-particle spectrum of the theory. At the same time the Keldysh propagator describes 
the state of the theory. Thus, our main concern below should be the solution of the DS equation 
for the Keldysh propagator . 

Due to spatial homogeneity of PP and due to its rapid expansion, which is supposed to fade 
away any initial inhomogeneity, we find it convenient to perform the Fourier transform of all 
quantities along the spatial directions: Dp'^'^{r]i,r]2) = J d^~^xe^^^D^'^'^{r]i,x;T]2,0). Then 
the Fourier transformed form of the DS equation for is as follows^: 



dr]3dri4 



Df(r/i,r/2) = <(r7i,772) + 
(27r)^-i J J {V3V4)° ^3) Df{r]3, r]^) D^_g{r]3, r]^) D^{r]4, 172) + 

+2 D^{r],, 773) Df{7^,, 774) ^4) (^4, m) + 2 /^f (r/i, I^f (^3, Va) D^-^iVs, V^) (^4, V2) + 

Note that here we use more symmetric form of the DS equation than the one in ([3]). In 
particular, under the integral, which defines loop contributions, we take both propagators on 
the sides of the self-energy E to be exact. 

We propose the following ansatz to solve ([5]): 



D {rii,ri2) = {r]ir]2) ^ d {pr]i,pr]2 



d {pr]i,pr]2) = -h{pr]i) h*{pr]2) 



1 + 2n{pT]i2) 



+ h{pT]i)h{pr]2) K[pr]i2) + c.c. 



(6) 



where 7712 = As well we use the tree-level retarded and ad- 

vanced propagators D^{r]i,T]2) = 6* (r/2 - r^i) (771772) ^ d' [pr]i,pr]2) , (7/1,7/2) = 

D-l , , , > 

-6^ (7/1 - 7/2) (7/17/2) 2 d~[pr]i,pr]2), where d' [pr]i,pr]2) = 2lm[h{pr]i)h*{pr]2)]. In ([6]) 
n{pr]) and K,{pr]) are unknown functions to be defined by the equations under derivation. 



Feynman rules can be found in [2o| . 
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This ansatz is inspired by the following observations. The retarded and advanced Green 
functions can be found as classical objects if the spectrum of quasi-particles is known. The 
ansatz for the Keldysh propagator follows from the interpretation of n{mj) as the particle 
density, {apttp), and of K,{pri) as the anomalous quantum average, {apQ-p), |8|. We assume that 
in the future infinity n and k are independent of the spatial coordinates. Furthermore, due to 
the symmetry of the PP under simultaneous rescalings of its coordinates, t] ^ Irj and x — )■ / x, 
we expect that in the future infinity n and k should be functions of the physical momentum, 
pry, only: np{ri) = n{pri) and Hp^rj) = n^pr]). 

It is known in condensed matter physics that non-vanishing k signals that one have chosen 
wrong harmonics to describe the quasi-particle spectrum. As well for constant k one can always 
set it to zero by performing Bogolyubov transformation which leads to the same ansatz ([6]), 
but with harmonics corresponding to a different a-vacuum and different value of n. We do not 
specify harmonics until the end where we check the IR behavior of K,{pr]) for the various choices 
of them. 

For general values of rji and ri2 the ansatz (jH]) does not solve the DS equation in question. 
However, in the limit ^771,2 — ^ and f]i/ri2 = const, one can neglect the difference between t]i 
and ?72 in the expressions which follow. That can be done if one keeps track only of the leading 
divergent IR contributions. 

As the result one can substitute the average conformal time 7712 = y/ViV2, instead of rji 
and r)2 separately, for the limits of integrations over 773 and 174. Then the ansatz in question 
reproduces itself under the substitution into the DS equation if n and k obey: 



(27r)^-i 



»?12 



dr]3 dr]i {r]3r]^) 2 x 



[ qr]3, qri4 ] d"^ [ \p - g|r/3, \p - g|r/4 ] +^d [ qrjs, qr]4\d [ \p - g|r/3, \p - g|r/4 



d [ qVs, W4 ]d^[\p- q\r]3, \p - q\Vi 



1 + 2n[pr]is) 



h*{pm) h[pr]i) + 



+4 e {t]4 -r]3) d {\p- q\T]3, \p - q\r]4 Re 



d [qvs, QVi ] h{pri3) h{pr]i) ^(^7742) 



(7) 
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and 



Vl2 



d [qV3,qV4]d [ \p - q\V3, \P - QlV-i ] + d [qV3,qV4]d [ \p - q\V3, \P ' Q\V4 ] + 



d [qV3,QV4]d [ \p - qlvsAP - llV-i 



1 + 2n{pr]ri) 



h* {pr]s) h* [prji) + 



+4 6'(?74 - 773) d^y\p-q\rj^, \p - q\riij d yqr]3,qVij h*{pr]3) h[prji) K.{prji2) | (8) 

where np°^ and n^p'' define the initial propagator (771, 772). Their presence is the drawback of 
the integral form of the equations under consideration, because then the equation itself depends 
on the initial conditions. The integrodifferential form of these equations is just the system of 
kinetic equations for n and k together, which was derived using different methods in |8|. 

In the derivation of ([7]) and ([8]) we have used the following relations d~{prii,pri2) = 



-d {pr]2,PVi) = - 
we assumed that n 



d~ {jyrji,prj2) and ^ d^~^qf {q,\p — q\) = f d^~^qf {\p — q\,q). As well 
prj) and ^(pr/) are slow functions in comparison with h{pri). Then one 
can safely change their positions under the drj^ and drj^ integrals, what we frequently do in 
the equations below. This is due to the usual separation of scales, which lays in the basis of 
the kinetic theory jl8|. In our case this approximation is correct only for the fields from the 
principle series, m > {D — l)/2, for which the harmonics h{pri) oscillate at the future infinity. 

However, the ansatz ([6]) solves ([5]) as well for the scalars from the complementary series, 
m < {D — 1) /2. For them the harmonics do not oscillate at future infinity. The main problem 
with the situation when h{pri) is as slow as n{pr]) and K,{pr]) is that then one can not derive the 
kinetic equation of the usual form. More complicated integrodifferential equations are available 
whose solution and physical interpretation is not yet known to us. 

To simplify ([7]) and ([8]) we change the variables as ^771,2,3,4 = 2^1,2,3,4 and use some approxi- 
mations [8j to arrive at: 
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n{pr] 



12 



n, 



(0) 



{2n 



iD-l 



1/v 



dq 



— / / dxs dX4, (X3 X4) 2 X 



d ix3,X4) 



+ 



n 2 



d {X3,X4) 



— d (X3, X4) d (X3, X4) 



p 



p 



X3 ] h { -X4 1+46' {x4 — xs) d (xs, X4) Re 



d (x3, X4) h ( -X3 ) /i 



1 + 2 n 



X4 K 



-Xl3 



P 



-X42 



and k(j9?7i2) 



,(0) _ 



Sd-2 



(2vr; 



1/7? 



12 dq 







/ / (ix3 dX4 (X3X4) 2 X 




-\ 2 



(0:3, X4) 



+ 



d (X3,X4) 



x/i* ( -X3 ] h* ( -X4 1 +4 6'(x4 — X3) (i^(x3,X4) (i (0:3, X4) h* ( -X3 1 h{ -X4 1 k 



+ (i {xs,X4) d (0:3, X4) 1 + 2 n 



p \ , I P 



P 



-Xl3 



1 



-X42 



Here is the volume of the {D — 2) -dimensional sphere of unit radius and Xij = ^Xi Xj . In 
(|9]) we have neglected p in comparison with q inside the integrals to keep only the leading IR 
terms. See js] for more detailed discussion. 

Now for the BD state h{x) oc 'Hf^{x). Then the ^3,4 integrals are saturated around x ~ /x, 
because of the rapid oscillations of the Hankel function at large values of their arguments. 
Hence, h{px3^4/q) can be expanded around zero, because p/q 1 in (Q. Then, because 
T-ll^Jix) behaves as C+ x*^ + C_ x~*'^, when x — )■ 0, there are interference terms under the dq/q 
integral which do not depend on q. As the result, both n{pr]) and K,{pr]) behave as A^ log(p?7) 
in the future infinity. Moreover, K{pr]) is generated even if it was set to zero at the initial stage 

k 

One should be a bit more careful with the similar manipulations for the other a-vacua, 
because their harmonics behave as linear combinations of e*^^ and e~*^^ at large momenta. But 
the careful study reveals the same picture for the most of them jsj. The explanation comes 
from the fact that their harmonics h{x) as well behave as linear combinations of x*'^ and x~*'^ 
in the future infinity. 

Only for the out-Jost harmonics, h{x) oc Ji^(x), which behave as single waves x*'^, the 
situation is different. In particular if one puts K{pi]) to be zero it is not generated back in 
(|9]). Or more precisely, contribution to it behaves as A^, i.e. is negligible in comparison with 
A^ log{pri). That is because the integrand of dq/q, defining k in (|9]), contains only g-dependent 
terms and, hence, the corresponding integral is convergent. At the same time for the harmonics 



.(9) 
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in question n{pri) has contributions of the order of X^log{pri). 

All in all, out-Jost harmonics represent the proper quasi-particle states in the future infinity. 
Which means that for out-Jost harmonics the ansatz ([6]) with K,{pr]) = does reproduce itself 
after the substitution into DS equation. That is possible if one neglects terms which are 
suppressed in comparison with powers of X^log{pri). This is the argument which favors the 
interpretation that independently of the initial state at the past infinity of PP the field theory 
state flows in the future infinity to the out-vacuum with some density of particles on top of it 
8|. At least if the initial state is just a mild excitation above an a-vacuum. To support such 
a conclusion we are going to show in a moment that for the out-Jost harmonics K,{pri) indeed 
flows to zero in the future infinity if it was not zero originally. 

The kinetic equation is obtained from the one for n{pri) in (|9]), with K,{pr]) set to zero, via 
application of the differential operator to its both sides: 



D-2 



dn{x) 

d\og{x) ~ ~2(27r)^-V 



g-3 

3^3 



dxs X ^ 



dxA X 



4 -^4 



X 



X ^Re[x-"'V{x3)xfV*{x4)] ^{1 + n{x)) nixsf - n{x) {1 + n{xs)y^ + 
+2Re[xfW{x3)x^'^W{x4)] ^(xs) (1 + ^(xg)) (1 + - (1 + ^(xg)) ^(xg) + 

+ Re [x^^F(x3)x4'''F*(x4)] [(l + n(x3))'(l + n(x)) - n(x3)'n(x)] | (10) 

Here x = p7]u, V{x) = [h' (x) - - . . . ] and W{x) = [\h {x)f - - • • 

where h{x) = Ji^J.{x) with J being the Bessel function. Dots in these expressions stand 

for a finite number of terms with higher powers of l/|x|. The presence of such contributions 
makes the collision integral well defined after the Taylor expansion of h{px/q) and can be 
explained by the behavior of the out-Jost harmonics in the limit x — )■ oo. All this is clarified 
in b|. 

n 

This is exactly the kinetic equation which was derived in [8[. If one have started with a 
small density perturbation over an a-vacuum, he can expect that n{pr]) is small in the future 
infinity. As is explained in js] in this case ffTOj) degenerates into a renormalization group type 
differential equation. The latter one can be solved with the result: 
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n{pri) 



r2 



D-2 



(2vr; 



D-1 



D-2 



(2 7r)^-V 



dyy' 
dyy 



C {priY +1 

Viy) 



where C is the integration constant, which depends on the initial conditions. 

This solution has stable point ^ ~ g-2vr^j ^ yu ^ 1, which approximately annihilates 
the collision integral in ( ITOjl . The stable point is reached when the production of particles is 
equilibrated by their decay l&j. In fact, from the collision integral ( ITOl) it should be clear that 
Fi defines the decay rate of the scalar particle into two, while r2 defines the particle production 
rate. Note that log(p?2) is decreasing as we approach the future infinity and n{pri) is the density 
per physical volume |8|. 

What is the most interesting fact, from the perspective of the discussion above, is that the 
stable point in question does not depend on A. Furthermore, it is not hard to see now that the 
stationary state of the kinetic equation ( fTOi) does as well annihilate, modulo subleading terms, 
the RHS of (H. 

The last thing which we have to check is the behavior of K,{pr]) for the out-Jost harmonics 
if it was initially non-zero. We as well assume that we have started from its small value in 
the past infinity and that it flows to the zero in the future. Under these assumptions, if one 
keeps only the leading terms the integrodifferential form of the equation for K,{pri) from 
degenerates to: 



d \og{pri) 



4:iX'^SD-2 ff\ , ^-if^ ^ 



{271 



dx^dx^x^'^ x^'^ """'^ 6'(x4 — X3) Im [V^(x3)V^*(x4)] . (12) 



00 



Here ReFs = Fi — F2 ~ (1 — e"^'"^) Fi > and, hence, the solution of this equation, K^prj) oc 
{P'lf^i flows to zero in the future infinity. Fe. our assumption is self consistent and ( fTTl) is 
stable under linearized perturbations of t^ijyrf). 
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IV. CONCLUSIONS AND ACKNOWLEDGMENTS 

We have found the result of IR dressing of a-vacua in PP of dS. The stationary state is 
described by out-Jost harmonics and corresponds to K{pr]) = with niprf) e~'^'^^. Further- 
more, the corresponding two-point correlation function depends, in the future infinity, only on 
the time difference, f]i/f]2 = e*^~*\ rather than on both of the times {rji and 772) independently. 

Using the same methods as those which lead to ffTOl) one can derive the kinetic equation in 
the contracting PP, ds"^ = dt^ — e~^* dx^ = ^ {drf — dx"^) where — )■ 77 = e* — )■ +00. The 
solution of the latter equation for low momenta p is jsj] : 



"^'^ ^ A-riog, ~ fi^' 

and is independent of p. It is valid for rj < r]Q = e™"'^*/'*'^ ^ 1. Here A is an integration constant, 
which depends on the initial state and f oc ^ > for m ^ {D — l)/2. 

One can see that the distribution in question grows with time, due to the contraction of the 
space and constant particle production, and moreover has a pole at some finite 770- In this case 
the backreaction on the gravitational background should be strong. This observation means 
that in global dS space the situation can be quite different form the one in the expanding PP, 
at least because global dS contains expanding and contracting PP simultaneously. Then we 
. lave two competing processes — expansion of the space time and explosive particle production 
8|. 
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RFBR-11-02-01120. This work was done under the support of the grant from the Ministry of 
Education and Science of the Russian Federation, contract No. 14.740.11.0081. 

V. APPENDIX 

The D-dimensional de Sitter (dS) space is the hyperboloid, = —Xq + Xf = 1, (/i = 
0,1, ... ,D and i = 1, . . . , D) in the {D + 1) -dimensional Minkowski space ds'^ = dX^ — dXf. 
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Throughout this paper we fix the curvature of the hyperboloid to be one. The expanding 
Poincare patch (PP) of this space is defined by the coordinates: 



Xq = sinh t + — e , Xd = — cosh t + — e' 

Xa = e'xa, a = 1, 1 (14) 

and covers only half of dS space, Xq — X^ = e* > 0. The induced metric in these coordinates 
is ds^ = dt^ — e^* dx"^ = ^ {diff — dx^), where rj = e~* = 1/{Xq — Xd). The past infinity of the 
PP corresponds to t — t- — oo, i.e. to ?7 = +oo. This is the boundary of the PP inside global dS 
space, Xq = Xo- The future infinity is at t = +oo, i.e. at ?7 = 0. The dS isometry is just the 
rotation symmetry group of the ambient Minkowski space, SO{D, 1). 

To quantize scalar fields in PP one has to specify time-dependent part of 

D-l 

the harmonics gp{i]) = ^'^2 ^ h(pr]) inside the harmonic expansion (/)(?7, x) = 
J d^~^p [dp gpi^i) e'"^^^ + fl'p(^) e*^^] . From the Klein-Gordon equation in PP follows that 



h{pri) has to solve Bessel equation with the index /i = ym^ — (■^^)^, where m is the mass of 
the particle. 

In time dependent backgrounds there is no basis of harmonics which can diagonalize the free 
Hamiltonian once and forever. The choice of the harmonics in the calculations corresponds to 
the choice of the background state \vac) = 0, usually referred to as vacuum. Bunch-Davies 
(BD) vacuum corresponds to h{pri) = e~ ^ T-lI^J {prj) , where H^^-* is the Hankel function of 
the first kind. These harmonics behave as e^^^ at the past infinity and diagonalize the free 
Hamiltonian only in that part of space-time. 

The other so called a-vacua can be obtained form the BD one via the corresponding Bo- 
golyubov transformations and correspond to the harmonics^ which are linear combinations of 
the Hankel functions of both kinds "H*-^^ and 'H^'^\ See e.g. |10] for a similar discussion in the 
global dS coordinates. 

Because of the time dependence of the Hamiltonian one has to apply the Schwinger-Keldysh 
diagrammatic technic instead of the Feynman one. In this technic every particle is characterized 
by the matrix of four propagators (see e.g. [l3], [isl): 
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G^(X,r) =z(</>(X)</)(y)), Gl_{X,Y)=i{<P{Y)<P{X)), 
G° +(X, Y) = (r0(X)0(F)) = e{T^y - G\iX, Y) + eiv, - r^,) G° _(X, F), 
G°_(X,F) = (f 0(X)0(F)) = ^^(r/, -r/.)G°_(X,F) + ^(r;,. -r^,)G%(X,F), (15) 

which obey one relation + G^_^_ = G"|^ + G^_. Here (T) T is the (anti-)time ordering. 
Note that the conformal time in our definition fiows in the reverse order (oo — ?■ r/ — )■ 0) with 
respect to the ordinary time t. 

All these propagators can be written with the use of the Whiteman function (^"(X, y) = 
i {(j){X)(f){Y)) . The latter solves the Klein-Gordon equation in the metric of PP. This equation 
is invariant under the full dS isometry although coordinates ( !T4|) are restricted only to the half 
of dS space. Hence, the solution of the Klein-Gordon equation should depend on the invariant 
distance between its two arguments — the two points on the hyperboloid, X^ = 1 and Y^ = 1. 
The convenient function of the latter one on the hyperboloid is the so called hyperbolic distance 
Z = -X„ F^. As follows from (HID it is equal to Z = 1 + in PP. 

The Klein-Gordon operator, when acting on the function of Z, rather than on the function 
of the two points X and Y separately, is equivalent to □((?) + = [Z"^ — l)i9| + D Z dz + to? 
|lo| . [q]. After the change of variables to x = (1 + Z)/2 the Klein-Gordon equation acquires 
the form of the hypergeometric one. Its solution (away from the singularity) is the following 
linear combination of the 2F1 hypergeometric functions: 



WW ^ . D-l D-l D 1 + Z\ 
G%Z) = A,F{ ^- + ^- - -; ^- ) + 



D-l . D-l . D l-Z\ / „ fD-V"^ 



+A2 F (^^— + — y ; J , ^i = \jm?-[ ) . (16) 

Here A12 are some constants which depend on the choice of the a-vacuum state with respect 
to which the averaging is done (see e.g. jlO|). For the BD vacuum A2 = 0. To take care of the 
behavior of this function at its poles and to obtain it as the quantum average (0(X)0(y)) one 
has to be more careful. 

In fact. Green function (fT6|) has three peculiar points in the complex Z-plane: Z = —1, 1, 00. 
They correspond to the usual peculiar points x = (1 + Z)/2 = 0, 1, 00 of the hypergeometric 
equation. The singular behavior G'^{Z) oc 1/{Z — l)-^/2-i corresponds to the situation when X 
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and Y sit on the same light-cone — the standard UV singularity of the propagator. Similar 
peculiarity of G^{Z) at Z = —1 corresponds to the situation when X sits on the light-cone with 
the apex at the antipodal point of Y. The antipodal point is obtained via the reflection at the 
origin of the ambient Minkwoski space jl^. Finally, at the infinity G^{Z) has the brunching 
point: lim^_j.oo G'^{Z) oc Z~^^ [Ci Z*^ + C2 Z~^^] with some constants Ci_2. 

To understand the behavior of G^{Z) at its poles it is instructive to consider Fourier trans- 
form of G^{Z) along the homogeneous spatial directions: 



-D-3 
2 



(0 {V.,P} -m ^ I d^'^'x e^^^'-y^ G\Z) = i!^:^ h{pv^) h*ipvy), (17) 



The appearance of different solutions of the Bessel equation in place of h{pri) here is in one- 
to-one correspondence with the concrete values of in flTB]) [lO|. 

Let us consider the BD propagator. Its only peculiarity inside the complex Z-plane is at 
Z = 1 and corresponds to the limit p — )■ 00 in momentum space. In this limit the Hankel 
functions behave as the plane waves. As it should be high momentum modes are not sensitive 
to the curvature of the space-time, i.e. they coincide with the flat space harmonics. For the 
inverse of the transformation ( |T71) to be well defined there should be an appropriate shift as 
Vx — Vy Vx — Vy ^ ^ ijT^ - The sign of this shift depends on which one among rjx and rjy 
is grater. As the result for the BD state (ill]: 



Gl4Z]=G'[Z-ze], Gl_[Z] = G'[Z - z e sgn{r]x - r^y)], 
G^_[Z] = G'^[Z + te], G^4Z] = G'^[Z + tesgn{r],-r]y)]. (18) 

Here G°(Z) is analytic on the complex Z-plane with the single cut going from Z = 1 to infinity 
along the real axis. 

The situation for the other a-vacua is different because in those situations harmonics are 
linear combinations of the Hankel functions of the two kinds T-L^^^ and 'H^'^\ The latter behave 
at large momenta as e"*^'' instead of e*^''. As the result for the other a-vacua (^^(Z) is defined 
on the complex Z-plane with two cuts connecting Z = 1 and Z = —1, correspondingly, to 
infinity and going, due to the ie shifts, in the opposite halfs of the complex Z-plane. 

Let us say a few words about the one-loop contribution to the propagators due to the A 0^ 
self-interaction. In the Schwinger-Keldysh diagrams there are two types of the vertices: of the 
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"+" and "-" type correspondingly. In the "+" ("-") type vertex only "+" ("-") ends of the 
propagators can terminate. Correspondingly the one loop correction can be written as: 

G\Zxy) = X' j [dW] j [dU] G\Zxw) G\Zuy) (19) 



where 



G°'i(Z)= ^ , and S°(Z)= L '\ L ^\ (20 

and the measure is [dW] = d^^'^^^W 5 (W^ — l) 9 {Wq — Wd), which is equivalent to the mea- 
sure ^ d^~^x on PP. This formula for G^ is valid for any a-vacuum. Note that in f ll9p dS 
isometry is naively broken by the presence of the Heavyside ^-function in the integration mea- 
sure, which restricts to the PP. 
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